The Circumbinary Outflow: A Protostellar Outflow Driven by a 

Circumbinary Disk 



Masahiro N. Machida 1 and Shu-ichiro Inutsuka 2 , and Tomoaki Matsumoto 3 

ABSTRACT 

The protostellar outflow is star's first cry at the moment of birth. The out- 
flows have indispensable role in the formation of single stars, because they carry 
off the excess angular momentum from the centre of the shrinking gas cloud, and 
permits further collapse to form a star. On the other hand, a significant frac- 
tion of stars is supposed to be born as binaries with circumbinary disk that are 
frequently observed. Here, we investigate the evolution of a magnetized rotating 
cloud using three-dimensional resistive MHD nested-grid code, and show that 
the outflow is driven by the circumbinary disk and has an important role even in 
the binary formation. After the adiabatic core formation in the collapsing cloud 
core, the magnetic flux is significantly removed from the centre of the cloud by 
the Ohmic dissipation. Since this removal makes the magnetic braking ineffec- 
tive, the adiabatic core continuously acquires the angular momentum to induce 
fragmentation and subsequent binary formation. The magnetic field accumulates 
in the circumbinary disk where the removal and accretion of magnetic field are 
balanced, and finally drives circumbinary outflow. This result explains the spec- 
tacular morphology of some specific young stellar objects such as L1551 IRS5. 
We can infer that most of the bipolar molecular outflows observed by low density 
tracers (i.e., CO) would correspond to circumbinary or circum-multiple outflows 
found in this report, since most of the young stellar objects are supposed to be 
binaries or multiples. 

Subject headings: binaries: general — ISM: clouds — ISM: jets and outflows — ISM: 
magnetic fields — MHD — stars: formation 
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Introduction 



Stars are born in molecular clouds which are magnetized interstellar clouds in our galaxy. 
Recent observations with high-angular-resolution measurements of polarized dust emission 
have shown that the st ar formation and ga s colla psing processes are certainly closely related 



to the magnetic field (IGirart et al. 1120061 . 12009T ) . The existence of highly c ollimated out 



flow/ j ets h a s established an ind ispensable role played by the magnetic field (ICabrit et al. 



19881 . 120071 ; iHirano et al. II2006I ). The realistic process of star formation in a magnetized 



molecular cloud is fully three-dimensional magneto-hydrodynamics (hereafter MHD) with 



merical simulation thanks to the recent development of various techniques ( 


kuc 


oh et al. 


1998; 


Tomisaka 


2002; 


Matsumoto & Tomisaka 


2004; 


Hennebelle & Fromane 


2008a 


). Hydro- 



dynamical calculations of protostellar collapse without magnetic field have shown that ini 



tially rotating molecular clouds tend to form binary or higher multiple stars (IGoodwin et al. 
20071 ). and the resultant binary stars are externally surrounded by circum-binary disks 



which have already been observed in star-forming regions. On the other hand, protostellar 
collapse calculations wit h magnetic field of the observationally inferred strength show spec- 
tacular bipolar outflows (lTomisakall2002l ; iMachida et al. II2004J ) . The outflow is driven by the 
Lorentz force and removes the angular momentum of the driving object, which suppresses th e 
fragmentation or the binary formation (IHennebelle fc Teyssierll2008l ; IMachida et al.ll2008bl ). 
Thus, the driving of MHD outflow and the binary formation have been thought to be in- 
compatible so far, although both the protostellar outflows and circumbinary disks have been 
observed numerously in star forming regions. Therefore, understanding of the puzzling link 
between them seems to require qualitative improvement in theoretical modeling. A clue for 
possible modifications are found in the evolution of magnetic field: most of the previous 
MHD calculations for the evolution of the collapsing cloud were done under the ideal MHD 
approximation that is valid only in the low-density phase of molecular gas but not in the 
high-density protostellar phase, which is obviously the room for the improvement for theo- 
reticians. Recently we studied the star formation process in resistive MHD calculations that 
account for the low ionization degree in molecular clouds, and showed the magnetic field is 
significantly removed from the centre of the collapsing cloud and frag mentation (or binary 
format ion) is actually possible even in a strongly magnetized cloud (IMachida et al.l 120071 
2008al lbL In this paper, we show the dramatic relation between the binary formation and 



outflow in resistive MHD calculations. 
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Model and Numerical Method 



To study the evolution of star-forming cores in a large dynamic range of spatial scale, a 
three-dimensional nested grid method is use d, in which the equat ions of resistive magneto- 
hydrodynamics are solved (see Eq. [l]-[5] of iMachida et al.l 120071 ). As the initial state, we 
take a spherical cloud with critical Bonnor-Ebert (BE) density profile, in which the uniform 
density is adopted outside the sphere (r > R c ). For BE density profile, we adopt the central 
density as n c = 10 6 cm -3 and isothermal temperature as T = 10 K. For these parameters, 
the critical BE radius is R c = 4.6 x 10 3 AU. To promote the contraction, we increase the 
density by 1.68 times. The mass inside r < R c is 0.8 M Q . The cloud rotates rigidly with 
Qo = 2.7 x 10~ 13 s _1 around the z-axis in the region of r < R c , while the uniform magnetic 
field (Bo = 32/iG) parallel to the z-axis (or rotation axis) is adopted in the whole computa- 
tional domain. In the region of r < R c , the ratio of thermal a^, rotational /?o, and magnetic 
7o to gravitational energies are a = 0.5, (3 = 0.04, and 7 = 0.06, respectively. 

For a realistic evolution of the magneti c field in the collapsing gas cloud, we adopted 
the resistivity (rf) derived in IMachida et al.l (120071 ). Figure [T]a shows the resistivity r\ and 
magnetic Reynolds number Re against the central number density (for detailed description 
see §2.2 of IMachida et al.ll2007l ). The gas temperature adopted in this study also plotted by 
the dotted line in this figure. 

To calculate a lar ge spatial scale, the nested grid method is adopted (for details see 
Machida et al. I l2005al ibl). Each level of a rectangular grid has the same number of cells 
(64 x 64 x 32). The calculation is first performed with five grid levels (I = 1—5). The 
box size of the coarsest grid I = 1 is chosen to be 2 5 i? c . Thus, grid of I = 1 has a box 
size of ~ 1.5 x 10 5 AU. A new finer grid is generated before the Jeans condition is violated. 
The maximum level of grids is / max = 13 that has a box size of 35 AU and the cell width of 
0.54 AU. 

We calculate two models: (a) resistive and (b) ideal MHD models. Hereafter, we call 
the former 'resistive model', and the latter 'ideal model.' Both models have the same initial 
condition shown above. The former include resistive term in induction equation, while the 
latter does not. 



3. Results 

As shown in Figure [T]a (the dotted line), the molecular gas obeys the isothermal equation 
of state with temperature of ~ 10 K until n c ~ 10 10 cm -3 , then cloud collapses almost 
adiabatically (10 10 cm~ 3 ^ n c ^ 10 16 cm -3 ; adiabatic phase) and quasi-static core (hereafter, 
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first core) forms during the adiabatic phase (jLarsonlll969l ; lMasunaga fc Inutsuka Il2000l ). The 



first core forms both in resistive and ideal MHD models when the central density reaches 
n c ~2x 10 12 cm -3 , and have a radius of ~ 10 AU and mass of 0.017 M & . By coincidence, 
the first core formation epoch almost agrees with the epoch at which the Ohmic dissipation 
becomes effective (i.e., Re < 1), as seen in Figure [T]a. Therefore, in resistive model, the 
magnetic field dissipates inside or around the first core with time after its formation. 

The evolutions of the magnetic field B zc and angular momentum J after the first core 
formation t c both in ideal and resistive models are plotted in Figure [T]6. The figure shows that 
B z _ c and J have the same evolution track in both models for t c ^ 100 yr but different tracks 
for t c ^ 100 yr. In ideal model, the magnetic field continues to increase to reach B z c ~ 1 G at 
t c ~ 1000 yr, while the angular momentum decreases to have J ~ 7 x 10 47 cm 2 s _1 at the same 
epoch. The decrease of the angular momentum is due to the magn etic braking which transfers 



the angular mo mentum outwardly along the magnetic field lines (IBasu fc Mouschoviaslll994 



Tomisakal 120021 ) . Since the cloud collapse slows and the rotation timescale becomes shorter 
than the collapsing timescale after the first core formation, magnetic field lines are strongly 
twisted and the angular momentum is considerably transferred owing of the amplified field. 
On the other hand, in resistive model, the magnetic field begins to decrease at t c ~ 300 yr, 
while the angular momentum continues to increase. Since the density of the first core exceeds 
n ^ 10 12 cm -3 , the magnetic field is effectively removed from the first core by the Ohmic 
dissipation as shown in Figure [T]a. Thus, the region around and inside the first core has 
a weak field that little contribute to the magnetic braking, thus, the first core in resistive 
model has a larger angular momentum than that in ideal model. As shown in Figure [U the 
magnetic field in resistive model is about two orders of magnitude weaker than that in ideal 
model, while the angular momentum in resistive model is about two orders of magnitude 
larger that that in ideal model. 

Figure [2] shows the distribution of density and magnetic field on z = and y = 
planes after the first core formation. The first core has a disk-like structure (Figs. [2]fc-n) 
with a centrally peak density profile (Figs. [2]a-rf). On the other hand, the magnetic field has 
a peak at edge or outside of the first core: the peak of the magnetic field has a ring-like 
structure as shown in Figures Wg-i- This is because the magnetic flux is removed from the 
region inside the first core by the Ohmic dissipation. Figures |5]&-g? indicate that the first 
core has a density of n = 10 12 — 10 15 cm -3 where the magnetic Reynolds number is below 
unity Re < 1 (see, Fig. [TJa), thus the magnetic field dissipates effectively. Since the first core 
increases in size with times, the decoupled region (i.e., the region with Re < 1) also expands. 
Therefore, the peak of the magnetic field strength moves outwardly as shown in Figures ^g-i. 
In resistive model, since the magnetic braking is not so effective owing to the weak field, 
the first core continues to spin up (see, Fig. []]&). In general, when the first core has a large 
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angular momentum, fragmentation occurs (e.g.. iGoodwin et al. 1120071 ). The fragmentation 
occurs ~ 800 yr after the first core formation (Fig. [2]e) only in resistive model. 

Both in ideal and resistive models, the outflow appears ~ 100 yr after the first core 
formation. This kind of outflow was already reported in man y past studies (e.g.. iTomisaka 
20021 ; iBanerjee fc Pudrital2006l ; iHennebelle fc Fromangll2008al ) . In Figures [2]fc- 0, the outflow- 
ing region is represented by the white-dotted line inside which the gas is outflowing from the 
center of the cloud. Outflow continues to driven near the first core in ideal model. On the 
other hand, in resistive model, the driving point of the outflow moves outwardly with time as 
seen in Figures |2]fc-o, because the decoupled region of Re < 1 expands outwardly. In resistive 
model, the magnetic field is too weak to drive outflow in the region near the first core where 
the magnetic dissipation is effective. Figure [3] left panel shows the configuration of the out- 
flow for resistive model, which indicates that the outflow driven by the circumbinary disk, 
not by the protostar (or proto-binary) . Figure [3] upper right panel shows the configuration 
of magnetic field lines that are strongly twisted by the rotation of the circumbinary disk. 
The circumbinary disk has a density of n < 10 12 cm~ 3 and is well coupled with the magnetic 
field that can drive the outflow. On the other hand, the magnetic field is too weak to drive 
outflow in the region around protobinary. Figure [3] shows B z on the z = plane with the 
same scale of left panel, and indicates that the magnetic field is weak around protobinary 
but strong in the circumbinary disk. 



4. Discussion and Summary 



In this study, we show that the outflow can be driven by the circumbinary disk. We 
calculated the evolution of outflow 1762 yr after its emergence, in which outflow continues 
to be driven by the circumbinary disk and reach 427 AU from the center of the cloud with 
a maximum speed of 8.4 km s -1 . At the end of the calculation, each fragment (i.e., each 
protobinary) has an only mass of ~ 0.03 M . During calculation, the protobinary revolves 
about the center of the cloud (i.e., common center of gravity) ~ 20 times, keeping a separation 
of about ~ 5 — 10 AU. It is expected that the mass of the protobinary, and size and speed 
of outflow increase with time in the gas accretion phase. 

Previous studies about the evolution of a rotating cloud showed that fragmentation fre- 
quently occurs only after the gas becomes optically thick (i.e., after the first core formation) 
and binary system is pos sible to form when the molecular cloud has a ce rtain amount of 
the angular momentum (IBodenheimer et al. 1 12000| ; IGoodwin et al. 1 120071 ). On the other 



hand, the study about a collapsing magnetized-rotating cloud showed that the magnetic 
field suppresses fragmentation and a single star tends to form in a strongly magnetized core, 
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because the rotation (or angular momentum) to promote fra gmentation is transferred by 
the magnetic effect such as the magnetic brakin g and outflow (IMachida et al. 112004 . l2005bl ; 
Price fc Batdl2007l ; Iffennebelle fc Teyssierll2008l ). However, these studies adopted ideal MHD 
approximation. In reality, the magn etic field dissipates in the high-density region in the col- 
lapsing cloud (INakano et al. II2002I ). Note that although the ambipolar diffusion and the 
Hall term effect can be also important fo r the magnetic dissip ation, the Ohmic dissipation 
dominates in the high-density gas region (IMachida et al.l 120071 ). 



We calculated the magnetized rotating cloud including the magnetic dissipation, and 
found the following: (1) fragmentation occurs to form the protobinary by the rapid rotation, 
since a weak magnetic field realized by the Ohmic dissipation makes the angular momentum 
transfer by the magnetic effect ineffective. In addition, such weak field cannot drive outflow 
from the protobinary, (2) strong outflow is driven by the circumbinary disk, because the 
magnetic field accumulates in the circumbinary disk where the removal and accretion of 
magnetic field are balanced, then the strong field is realized. The binary system often shows 
a single outflow that may be driven by the circumbinary disk, not by one protostar in binary 
system. For example, a single molecular outflow in L1551 IRS5, in which binary system is 
embedded, may be driven by this mechanism. 

In the present study, we adopted the resistivity rj as most plausible value as shown in 
Figure [T]a. However, the resistivity depends on the size distribution and abundance of dusts. 
The resistivity is closely related to the magnetic dissipation. The magnetic field is related to 
most fundamental phenomena of the star formation: the angular momentum and magnetic 
flux problem, fragmentation and binary formation, and outflow driving. Thus, it is crucial to 
determine the resistivity for resolving major problems in star formation. Since the strength 
and driving point of the outflow is directly related to the resistivity as shown in this paper, 
we can determine it with future higher-resolution observations such as ALMA. 



Numerical computation in this work was carried out at the Yukawa Institute Computer 
Facility. This work was supported by the Grant-in- Aid for the Global COE Program "The 
Next Generation of Physics, Spun from Universality and Emergence" from the Ministry 
of Education, Culture, Sports, Science and Technology (MEXT) of Japan, and partially 
supported by the Grants-in-Aid from MEXT (18740104, 20540238, 21740136). 



REFERENCES 

Cabrit, S., Goldsmith, P. F., & Snell, R. L. 1988, ApJ, 334, 196 



-7- 

Cabrit, S., Codella, C, Gueth, F., Nisini, B., Gusdorf, A., Dougados, C, & Bacciotti, F. 
2007, A&A, 468, L29 

Girart, J. M., Rao, R., & Marrone, D. P. 2006, Science, 313, 812 

Girart, J. M., Beltran, M. T., Zhang, Q., Rao, R., & Estalella, R. 2009, Science, 324, 1408 
Banerjee, R., & Pudritz, R. E. 2006, ApJ, 641, 949 
Basu, S., & Mouschovias, T. C. 1994, ApJ, 432, 720 

Bodenheimer P., Burkert A., Klein R. L, k Boss A. P., 2000, in Mannings V., Boss A. P., 
Russell S. S., eds, Protostars and Planets IV. Univ. Arizona Press, , p. 675 

Goodwin S. P., Kroupa P., Goodman A., k Burkert A., 2007, in Reipurth B., Jewitt D., 
Keil K., eds, Protostars and Planets V. Univ. Arizona Press, , p. 133 

Hirano, N., Liu, S.-Y., Shang, H., Ho, P. T. P., Huang, H.-C, Kuan, Y.-J., McCaughrean, 
M. J., k Zhang, Q. 2006, ApJ, 636, L141 

Hennebelle, P., k Fromang, S. 2008a, A&A, 477, 9 

Hennebelle, P., k Teyssier, R. 2008b, A&A, 477, 25 

Kudoh, T., Matsumoto, R., & Shibata, K. 1998, ApJ, 508, 186 

Larson, R. B., 1969, MNRAS, 145, 271. 

Machida, M. N., Tomisaka, K., & Matsumoto, T. 2004, MNRAS, 348, LI 

Machida, M. N., Matsumoto, T., Tomisaka, K., & Hanawa, T. 2005a, MNRAS, 362, 369 

Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2005b, MNRAS, 362, 382 

Machida, M. N., Inutsuka, S.-L, & Matsumoto, T. 2007, ApJ, 670, 1198 

Machida, M. N., Inutsuka, S.-L, & Matsumoto, T. 2008a, ApJ, 676, 1088 

Machida, M. N., Tomisaka, K., Matsumoto, T., & Inutsuka, S.-i. 2008b, ApJ, 677, 327 

Masunaga, H., & Inutsuka, S., 2000, ApJ, 531, 350 

Matsumoto, T., & Tomisaka, K. 2004, ApJ, 616, 266 

Nakano, T., Nishi, R., & Umebayashi, T. 2002, ApJ, 573, 199 



- 8- 

Price, D. J., & Bate, M. R. 2007, MNRAS, 377, 77 
Tomisaka K. 2002, ApJ, 575, 306 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 9- 




Fig. 1.— (a) Resistivity (77; left axis), magnetic Reynolds number (Re; right axis), and 
temperature (T; outer right axis) as a function of the number density. The magnetic field 
well couple with the gas in "coupled" region, while magnetic field decuples from the gas 
in "decoupled" region (grey zone). The first core formation epoch is plotted by the arrow, 
(b) The magnetic field strength B z c (left axis, black lines) at the center of the cloud, and 
angular momentum J (right axis, red lines) in the region of n > 0.1 n c against the elapsed 
time after the first core formation for ideal (dotted lines) and resistive (solid lines) models. 
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Fig. 2. — Distributions of density (panels a-e, and k-o) and magnetic field strength (panels 
f-j and p-t) around the first core (grid level / = 13) on z = (panels a-j) and y = (panels 
k-t) plane. Arrows in panels a-j mean the velocity vector, while those in panels k-t mean 
the magnetic field vector. The elapsed time from the initial (t) and that (t c ) after the first 
core formation, and maximum number density (n max ) are described in each top panel. The 
dotted line in panels k-o means the border between infall and outflow inside which the gas 
is outflowing from the center of the cloud. 
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Fig. 3. — Large-scale structure of I = 10 grid in three dimension. Left panel: The config- 
uration of circumbinary outflow is shown by red-volume, in which color means the outflow 
speed. The density contours (colors) and velocity vectors (arrows) are projected in each 
wall surface. The square at lower right corner shows the close-up view around the proto- 
binary, in which color means the density distribution, and the high-density region (i.e., the 
protobinary) is represented by the red iso-density surface. Right upper panel: The magnetic 
field lines integrated from the circumbinary disk are plotted. The bottom panel indicates 
the strength of B z by the color. Right lower panel: The strength of B z on z = plane is 
represented by the color. 



